


<n. ..?> 




rj;k% isc.,'; 
£.r-..:v- 


54^5 X^,.', IFSH- 73 -Iq 


; v<€^X ’■ 


LEHIGH UNIVERSITY 


ll/ffTTllkl 


\ A /x* 


•'./■ "X 




.ms 


IftP# ISp BONDED HALF’ PLANES CONTAIN I NG •■•;/ - h • . 

Jr-./-. - ; : .-.v. .{■. .JJ .' 

AW^ADDTTDADTr V HOT ETMTCTl TDAPI/ '- t&WVs. 


mmm 



V W-iawxM 


./V. !<..■*■■■%?:>:■< ' v-fcS^i 
- is it •' • s .'. n- 


:Y ■’Sr. ' X Tl-n^lX •; -.'■>■ X- X ~<Sl X* v i ;; 'V \ - - 

-.'Av -z. -•• ! :•• : ., • 1 » * •..". •; • .. . >.\U : ■ . r.-iiV-X A-. 


r --.iv' ■ -/ • . / •■•■^>V‘ • ■,44 i,^WvA- 

XXv ; ./\ i - j ';'-J:. : ": ;.. ".r V.. H' . -. . : ;. .. :' : 


; f V ' Technical 1 Report NASA TR- 7 J -9 

* ‘ k' * ' . «. • • •■-/.. 


'••••■.m- vv;m: 


• w**'^.* .-*-4 V V '‘ N - 

: 1 )i • . t *' *” _« ft. U, 


:•)} ; . > 

ViilVl te . 


,xy- 

; !ll* . , i L - !).C 


..’ • •' ‘ -' > / •••’■V.T ^ ■' 

. vwf‘,-' / < :; - .« 


-V 'I 


■>■■■•■ ■ V 


isvr 

• ' i 




; - 

x* 'fi a V- 


,r>i H 


r^r::-- 7 






>-• 

1 




m ' si- ^ 

Its ' . '*•£%&, 9 *> s 




\\ 

I t[r hi H 

\ ( >>... !i ‘n 
K ’ 


Iay 1973 


ft 


a National Aeronautics and Space Administration S 4 - 'x 


‘•A ■ : r.\’ i- '. : . ' T\~^ 7 fy*-l:y -. 

•Grant NGR-39-007-011 „ • ' 

J -. ........ f .:•■■ . 


: — ” T s. 

«•- / < - v/.i? ^ 


... 

feuT ‘ 1 la 


h£S-' : n- v - V- : --r-vV/r-:- 

y.- • : 


-f" ';/ 


-*M.ri-KrS • 


•r.»t . -■ • • -|.*v V; 

' . .-V •*V», ;.\i > . -.• :■ f •• -.^'.rv- •]■ ; ■ i/v: , . 

. .',^v '.?■ . ' \; ■ v ■ ■/; -V • jr ;.. r - 


-‘-rVKf .r *V 

: 

r-X' ** • 


mssv ; x3Mxi\ : . 

*//,/ . . • . -; 4 » fit : ^ >: 

■ T. - rr if. \ - 


. *v % - . 4 * , . rn v- - * * :> 4 . 

c ■• - 1\\' jx " !'* ’ 


■ \ 


V V 1 - . 




BONDED HALF PLANES CONTAINING 
AN ARBITRARILY ORIENTED CRACK* 

by 

F. Erdogan and 0. Aksogan 
Lehigh University, Bethlehem, Pa. 


Abstract 

The plane elastostatic problem for two bonded half planes 
containing an arbitrarily oriented crack in the neighborhood of 
the interface is considered. Using Mellin Transforms, the prob- 
lem is formulated as a system of singular integral equations. 

The equations are solved for various crack orientations, material 
combinations, and external loads. The numerical results given in 
the paper include the stress intensity factors, the strain energy 
release rates, and the probable cleavage angles giving the direc- 
tion of crack propagation. 


1 ■ INTRODUCTION 

The structural strength of composite materials is controlled 
to a considerable extent by the size, shape, orientation, and 
distribution of the flaws and imperfections which exist in the 
material. Usually these flaws and imperfections exhibit them- 
selves in the form of entrapped gas or weak impurities on the 
interface, ruptured bonds, cracks, inclusions, and geometric 
s i ngul ari ties arising from the particular shape of the constitu- 
ent materials. From the viewpoint of fracture initiation and 
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propagation in the medium particularly important are the manufac- 
turing flaws such as flat cavities which develop during bonding 
or casting, small cracks resulting from the residual stresses, 
and fatigue cracks caused by the cyclic nature of the external 
loads. Thus, in studies relating to the fracture initiation and 
propagation in the material, it is necessary to have a good esti- 
mate of those factors representing the severity of the external ' 
loads in the neighborhood of the "isolated dominant flaw". Gener- 
ally, the comparison of these factors (such as the stress inten- 
sity factors, the strain energy release rate, the crack opening 
stretch, or the cleavage stress at a characteristic distance from 
the flaw boundary) with the corresponding characteristic constant 
representing the resistance of the material to fracture consti- 
tutes the fracture criterion. If the dominant imperfection is 

completely imbedded in a homogeneous phase and is located suffi- 
ciently far from the phase boundaries or interfaces, then the 
disturbance of the stress field around the imperfection will not 
be affected by the neighboring phases and the disturbed stress 
field may be obtained by solving the problem for an infinite 
homogeneous solid. On the other hand if the flaw is located near 
a phase boundary or a bi-material interface, then the solution 
of the problem for the nonhomogeneous medium becomes necessary. 

An up-to-date review of the available solutions for variety 
of crack and inclusion geometries in composite materials may be 
found in [1]. The primary interest of this paper is in the 
evaluation of the disturbed stress field around a crack located 



near a bi-material interface or, as a limiting special case, 
free boundary. The problem was studied in a previous series of 
papers for two special crack orientations, namely- the case of a 
crack located parallel to or at the interface [2-6] and the 
problem of a crack perpendicular to and crossing the interface 
[7, 8]. In this paper we will assume that the orientation of 
the crack with respect to the interface (i.e., the angle 0 Q and 
the distance d in Figure 1) is arbitrary and the problem is one 
of plane strain or generalized plane stress. As in the previous 
studies it will also be assumed that the interface in the non- 
homogeneous medium is either a plane or has a sufficiently large 

' i - ' 

radius of curvature so that the disturbed stress field can be 
approximated by that of a crack in two bonded elastic half 
planes (Figure 1). Even though the problem will be formulated 
and solved for the bonded half planes, the technique described 
in the paper appears to be quite general and may be used to treat 
the problem of any number of bonded wedges with radial cracks. 

2 . THE INTEGRAL EQUATIONS OF THE PROBLEM 

Using the conventional superposition technique the solution 
of the problem of a traction-free crack in the composite medium 
under a given set of external loads can be expressed as the sum 
of two solutions: the first obtained for the given external loads 
and the given medium without the crack, and the second obtained 
for the two bonded half planes with a crack where the only exter- 
nal loads are the crack surface tractions which are equal and 
opposite to the stresses found in the first solution on the 




presumed location of the crack. It is clear that only the second 
solution which gives the disturbed stress field due to the exist- 
ence of the crack will have singularities. Also note that in the 
second problem, since the external loads are local and statically 
self-equilibrating, in the application of Mellin transforms the 
regularity conditions required of the solution as r -»■<»' will be 
satisfied. Following now, for example [9], in polar coordinates 
the plane elasticity problem for a medium having the elastic 
constants y and k (k = 3-4v for plane strain, k = (3-v)/(l+v) for 
plane stress) may be formulated as 


a(r,0) = x rQ + i t Q Q 
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Referring to Figure 1, the medium will be considered as 
consisting of three infinite wedges: 1 (y-| >k-| ) , < 8 < *> 


2 ( ^ 2 ’ K 2 ^ ’ e 0 < 9 < 2 * ^ ~ 

wedge with an arbitrary angle, using the Mellin transforms to 
solve (2), from (1) we obtain 


•nr 


2 <9<6 0 


For an infinite 
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M[r 0 ,s] = E(s,6) 


= 2i(s+l)[Ase 1se '+ Bfs + De 1 *^ 2 ' 0 - Be- 1 ' 5 * 2 ’ 6 ] , 


M[r 2 T rr> s] = - s(s+l)(Ae is0 + Ae‘ ise ) 


- (s + l)(s + 4)[Be i(s+2 > 6 + ffe- 1(st2)0 ] , 


M[r^v,s] = V(s,6) 

[Ase isn + B(s-M)c i ( s+2)B + kSc-''( s+2 ) 0 ] , 

/ (3 .a-c) 

where A and B are functions of the transform variable s. The 
Mellin transform of a function f(r) defined and suitably regular 
in (0<r<«>) s and its inverse are defined by 

F ( s') = M [ f j s ] = / f(r)r s_1 dr , . - 


f(r) = / F( s ) r -:> ds , 

c-i°° 


c+i 00 


-s. 


{A. a , b ) 


where c is such that r c_ ^f(r) is absolutely integrable in (O, 00 ). 
The transform of derivatives may be shown to be 

0 dr' 
provided 


J r n d^lrl r s-l dr = ( . 1} n Ii|^l F (s) , 


(5) 


.s+m-1 d m ~ 1 f 


dr 


m 


.-*■ 0 as r ->• (0,“) , (m=l , . . . ,n) 


( 6 ) 


In applications conditions (6) provide the information to deter- 
mine the strip of regularity containing the line Re(s) = c in 
the inversion integral. 


It should be noted that in working with Mellin transforms up to 
the inversion^stage in the manipulations the transform variable 



If we now let the subscripts 1,2,3 stand for the wedges 
shown in Figure 1, the problem must be solved under the following 
boundary conditions: 

o-j (r,u/2) = o 2 (r,TT/2) , (0<r«») , 

v-| ( r , ir/2 ) = v 2 (r, 7 r/ 2 ) , (0<r<“) , (5.a,b) 

a 1 (r,3ir/2) = 0 3 (r,- 7 r/ 2 ) , (0<r<«>), 

v-j (r,3ir/2) = v 3 (r,-ir/2) , (0<r«»), ( 6 .a,b) 

cT 2 (r, 0 o ) = a 3 ( r > 6 0 ) > (0<r<«>), (7) 

lim cr 9 (r,e) = p 9 (r) + ip-,(r) , (a<r<b), (8) 

0-0 o+ O c c 

v 2 (r,0 Q ) = v 3 (r ,0 o ) , (0<r<a , b<r«=°) , (9) 

/ C v 2 ( r ’ e 0 ) ~ v 3 ^ r, 0 o^ dr = 0 * 0°) 

where 0 j and,y . , (j = l,2,3) are defined by (1), y 2 =y 3 , and note 
that (9) and (10) correspond to the condition of displacement 
continuity outside the crack. In the solution of the problem as 
given by (3) there are six unknown (complex) functions A-(s), 

B • ( s ) , ( j = 1 ,2,3) to be determined. The homogeneous conditions 

J 

( 5 ) - ( 7 ) provide five equations. The sixth equation is obtained 
from the mixed conditions given by ( 8 ) and (9). Thus, eliminat- 
ing five of the unknown functions, the problem may be formulated 

s is treated as a real variable. The function F(s) is analyti- 
cally extended into the complex plane from the real line when 
the inversion is evaluated. In the present problem, the complex 
notation in (1) and (3) is used only for convenience. Thus, for 
example, in separating the transforms of the stress components 
T rg and tgg in (3. a), s should be treated as real and A and B 

should be treated as complex quantities. 
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as a system of dual integral equations for tv/o unknown real 
functions by using (8) and (9). However, a somewhat more direct 
method to solve the problem would be its reduction to a system 
of singular integral equations for a pair of real’ functions f-j 

■k 

and f 2 defined by 

v 2 (r,0 o +O) - v 3 ( r > 9 0 - 0 ) = f(r) = f 2 (r) + if-,(r) , 

(a<r<b) . 

If' (8) and (9) are replaced by 

ff(r) , (a<r<b) . , 

v 2 (r,e 0 ) - v 3 (r,e 0 ) 
and i f we defi ne 


0 , (0<r<a, b<r<«>). 


( 11 ) 


( 12 ) 


U(s) = M[r 2 {v 2 (r,0 o ) - v 3 (r,6 0 )},s] = J f(r)r s+1 dr 

a 


u c + l 

= J [f 2 (r) + if l( r )] r dr = u 2 (s) + iu l( s ) 

3 


(13) 


by substituting from (3) into ( 5 ) - ( 7 ) and (13) we obtain six 
linear algebraic equations in A i and B - , (j = l,2,3), which may be 
solved giving A-(s) and B-(s) in terms of U(s). For example, 

J J 

for the wedge 2 we find 


Ao(s) = 


- e 


(1+k 2 )s(s+1 ) (e 
+ m,(s + l) 2 e 1(s+2 > e o 


nrs „-vrrs 


{[m-|(s + l)se ls0 ° 


_ se- is(7r+0o) + m 2 e i(s+2)6 °]U 2 (s) 


+ [m 1 (s+l)(s+2)e is0 ° + m ] (s+1 ) 2 e i (s+2)0 ° 

+ (s+2)e" is(Tr+8o) + m 2 e i(s+2)0 °]iU 1 (s)} , 


Mote that physically f] and f2 represent the densities of edge 
dislocations distributed along 0 = 9 0 in two bonded half planes. 
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B 2 ( s ) = 


( 1 +k 2 ) (s + 1 ) (e" 17rs - e 17TS ) 

- mi (s + l)e i(s+2)0 o - e ~T( s0 o + 29 0 + s ' n ')] U ^ (s) 

- [m 1 (s+2)e is0 ° + m 1 (s+l)e i(s+2)0 o 
1 (se 0 + 20 o + STT) ]fu ^ (s)} } 


{ [-m-j se 


i s © o 


- e 


(14. a, b) 


where 


m-j = (m-1 )/ (ic^m+l ) , 

m 2 = ( 1 <-|-mK:2}/(»<i +m ) 


m ~ y -j / y 2 * 


(15) 


Expressions similar to (14) may be found for wedges 1 and 3. 

All the field quantities in wedge 2 may now be obtained in 
terms of f(r) by substituting from (14) into (3) and using the 
inversion formula. In particular, from (14), (3. a), (8), and 
(4 . b) we obtain 

1 +k 9 , b 

Cp 2 (r) +.1p,(r)]- lim 1 / tK 2 (r.r 0 .e)f 2 (r 0 ) 


9+9 o +0 a 


+ iK i ( r » r 0 > e ) f i ( r o^ dr o v (a<r<b) , (16) 


K 2 (r>r 0 ,9) = 


c+i°° . r s+1 

f ds r o 

J . ins -ins s+2 
c-i°°e -e r 


{ro,s(s+l}e 1s(9+e ° ) 


+ [ro 2 + m ) (s+l) 2 ]e 1 ’ (se + s9o + 26o) . se f s < 6 ' e ° ' 1r > 

+" m,s(s + l)e its8 + s9o + 26) + mi(s+l) 2 e i(s+2)(9+9o) 
+ ( s + i ) e i (s 9 + 29 - S 0 O - 29 0 - stt) 

+ m ,se - 1(s9 + 2e+s9o) + m,(s-H)e- 1 ( 5 + 2 )( 9 +e o> 

+ e i (s 9 o + 29 0 - s 9 - 29 + sir) i 


K 1 ( r,r o’ 0 ) = 


c+i 00 

I 


ds 


s+1 


„ . ^i^s - i tt s s + 2 

c- 1 co e - e r 


?T5- {m 1 (s-M)(s + 2)e 1s(0+e '>> 


+ [m 2 + n. 1 (s+1) 2 ]6 1(se ° > + s 0 + 20o) 

' + (s+2)e 1 ‘ s(e - 9 °- lr) + m l( s+l)(s+2)e 1(se+ s8 ° +26) 

_ ( s+1 ) e 1(s 0 + 26 - S0 O - 200 - str) 

- m 1 (s+2)e" ,(s? + s8 ° + 26) - m 1 (s+l)e‘ 1(s+2)(e+eo) 
; + ;,1(s9 0 + 29 0 - S9 - 26 + sir) j / ( 1 7 - a , b ) • 


In (16) the order of integrations has been changed. For 0> 0 O , 
since the related integrals are uniformly convergent, this is 
permissible. (16) provide a system of integral equations to 
determine the unknown functions f-j and f 2 . To solve this system 
the kernels K-j and K 2 must be evaluated which may be done either 
by using the residue theory and expressing them as infinite series 
or by reducing the integrals to real integrals and evaluating them 
numerically. In either case, it is first necessary to determine 
the strip of regularity containing the constant c. Let the inte- 
grands in (17) be analytically extended into the entire plane and 
let the poles s. be ordered as 

J . 

--- < Re(s_ 2 ) < Re( s_i ) < Re(s +] ) < Re(s +2 ) < --- (18) 

Let Re(s_-j) < c < Re(s + -j). Thus, the integrals in (17) may be 
evaluated by closing the contour to the left for r < r Q , and to 
the right for r-> r by means of semicircles of infinite radius, 
and by summing the residues. Noting that for r->*0, T..~0(r"^), 

* vl 

Re ( X) <1, and for r t . . ~ 0 (r~ a ) , Re (a) >_ 1 , and since the 
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residues are of the form r~^ s J + ^, s. then should be ordered 
such that 

Re ( s_-| ) < -1 , Re(s +1 ) > -1 . (19) 

Thus, Re(s_^) < Re(s) = c < -1 gives the strip of regularity, 
s i being the first pole to the left of the line Re(s) = -1. 

For an arbitrary value of 0 , even though the residues in 
(17) may be evaluated without any difficulty, the resulting infi- 
nite series cannot be summed in closed form. Hence, it is diffi- 
cult to study the singular behavior of the kernels. For this 

i ’ . . , ■ ’ •• ' 

reason in this paper the kernels will be evaluated by reducing 
(17) to real integrals. For this we let c = -1 , s = -1 + iy, 
(~oo<y<oo), and indent the contour in such a way that the pole 
s + i - -1 lies to the right of the line of integration. The inte- 
grals in (17) may then be expressed as 

c+i co — ■ 

K k^ r,r o 5 °^ = / H k (s)ds 

C-i°° • ; . ■*“ V 

oo ; 

= / H.(-l+iy)idy - lim tri ( s+1 ) H. ( s ) , 

-00 S->-l • 

( k = l ,2). (20) 

Evaluating the residue at s = -1, from (17) it may be shown that 

lim Tri (s + 1 )H. (s) = (m ? - m, - 2) , (k=l ,2) . (21 ) 

s->-l K r c 1 

On the other hand, from (10) and (11) we have 

b 

/ f k (r 0 ^ dr 0 = 0 ’ (k=l ,2). (22) 

di . 

Thus, when substituted into (16), the integrated terms in (20) 
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will have no contribution. In the remaining integral in (20) 
the integrand turns out to be the sum of an odd function and 
an even function* in y, giving K. in terms of real integrals in 

J 

As r -»• r Q these integrals become divergent. Since the 
integrands are bounded at y=0, the divergent parts can be sepa- 
rated by considering the asymptotic behavior of the integrands as 
y -> «>. By defining 


P = log(r Q /r) , e = 6 - 0 Q , 


(23) 


for small values of e from (17) and (20) it may be shown that 


K k (r.r # .«) = 1; {e-^D + O( e -W- 20 °))]si n py 

+ 0(o” y ^ : '’ 2 ' 0+,: ^ JcospySdy 

* f + M k (r,r 0 ,e) , (O<0 o <J, k = l ,2) , (24) 

P +e 


where M k ( r , r Q , 0 ) 
closed interval 


is bounded for all 
[a 5 b] . 


values of r and r Q i n the 


Substituting now from (24) into (16), separating the real 
and imaginary parts, and letting e->0 we obtain 


1 +K, 

2m 0 


P-j ( r j 


= 1 J 
u a 


2 

I 

j = l 


6 . . 

1 J 


1 og(r Q /r) 


k ij< r >o )]f j^ )dr o 


( i = 1 ,2; a<r<b ) , (25) 

where the bounded kernels k.., (i,j=l,2) are given by 

' J 

00 

k ll< r * r o) = h l iTsfe s1npy[cosh2e o y(n, l + 4n, 1 J,2cos 8 o - m 2 ) 

+ 2m-jy sin20 o sinh20 Q y + 2e~ 7ry ] , 

" " 

If not, it can always be put in that form. 
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k 12^ r>r o^ 


J_ 7 [cospy sinh26 o y(4m 1 y 2 cos Z 0 o 

2r J Q si nhtry . 0 

. n - m 2 ) - 4 mi y cos2e o sinM s 1nh . z0 o y] ’ 


k 21 (r,r 0 ) = 


1 r Al — [cospy sinh20 y(m-| - 4m-jy cos 0 Q 
2 r J Q sinhivy L 0 


+ m 2 ) -■ 4 m^y cos 2 0 q sinpy sinh2e Q y] , 


, , - ; . l_-7 sinpy[cosh29 y (m-j + 4m 1 y 2 cos 0 Q - m 2 ) 

k 22^ r,r o } " 2r J sTnhiy . 0 1 


- 2m- 


sin28 rt si nh20 y + 2e * y 3 » 


(a < (r ,r ) < b) . 


(26 .a-d) 


In the system of integral equations (25) the dominant kernels 

have a simple Cauchy-type' singularity. This may be seen by 


observing that 


FTogTr 0 /rT- r iog[l + (—--1)3 


_J C , + jizlf (>-l) B l 

r(~°- - 1 ) 1' " +1 


n -1 


r o- r 


[l + o(^~- - D3 • 


(27) 


Thus, the system of singular integral equations is of the 
following conventional form: 


^Pi(r) riJ ! + h ij (r,r o nf j (r 0 )dr 0 . 

2u 2 1 Tr a j = l o 

(i = l ,2; a<r<b) , 


h ij(T,r 0 ) » kijtn.ro) + [rW(r 0 /ry ' r 0 -r ]6 ij ’ 


(1,0 = 1.2)- (28. a, b) 






The index of the integral equations (28) is ic = 1 ; hence the 
solution will contain two arbitrary (real) constants which are 
determined from the additional conditions (22). 

3. SOLUTION OF THE INTEGRAL EQUATIONS AND STRESS INTENSITY 
FACTORS 



Referring to [10], it may be shown that the index of the 
system of singular integral equations is +1 and its solution is 
of the following form: 

fj (r ) = gj (r) [ ( b-r ) (r-a)]" 1/2 , (a<r<b, i=l,2), (29) 

where the unknown functions g^ ( r) , (i=l,2) are bounded in the 
closed interval [a,b]. Even though, in principle, the system 
can be regularized and reduced to a pair of Fredhol m-type inte- 
graT equations , its solution may be obtained with much less com- 
putational effort by using the technique described in [11]. In 
the numerical solution the main problem is the evaluation of the 
kernels k.. given by (26) for which in this paper a modified 

* J 

version of Filon's integration formula has been used [12]. 

In the application of the results to fracture problems in 
composites, of particular interest. are the stress intensity 
factors and the probable plane of cleavage at a given crack tip. 
The normal and shear components, k-j and k 2 , of the stress inten- 
sity factor are defined by and may be evaluated from the follow- 
ing expressions [2, 4, 5]: ' 

2p- ' 

k i ( a ) = lim /2(a-r) T 200 (r,6 o ) = lim yyy— /2(r-a) fj(r) , 






: 2 ]i ? 

k -j ( fa ) = 1 i m /2(r-b) T 2QQ (r-,e o ) = - lim /2(b-r) f^r) , 

■ 2y ? 

k ? (a) = lim /2(a-r) t« fi (r,0 )■ = Tim - /2(r-a) f ? (r) , 
d r->a ^ r0 0 r+a 1 K 2 / 

. 2y„ 

k ? (b) = lim /2(r-b) t ? fi (r,6 ) = - lim yr— ■ - V2(b-r) f ? (r) . 

6 r+b ^ r0 0 >>b - 2 . 

(30. a-d ) 

The constants k-j and k 2 are a measure of the intensity of 
the stresses around the crack tips. For example, at the crack 
tip (r=b, 0=0 O ) the cleavage stress may be expressed as [13] . 

^dxf) ( P a = cos 2 ( k i cos2 f - f k 2 s i n<J> ) + 0(/p) , (31) 

/2p 

where (p,<J>) are the polar coordinates at the crack tip, <j> being 
measured from the line which is the prolongation of the crack 
(Figure 1). Thus, once k-j and k 2 are determined, for brittle 
solids the probable angle'. 'of crack propagation may be postu- 
lated as the angle of the radial plane corresponding to the 
maximum cleavage stress and may be determined from = 0, 

? o 

and d o^^/dtp < 0, or 

k 2 ( 1 - 3cos<f> c ) - k-j si n<f> c = 0 , 

3k 2 sin<f) c - k -j c o s 4> c < O'. (32. a, b) 

If an energy balance type criterion is used to estimate the crack 
propagation load, one may need to calculate the strain energy 
release rate which is given in terms of k-j and k 2 as follows [13]: 

■ j( 4p K) < k i 2 - 1 - ■ (a 0 = (b-a)/2). (33) 


4. THE RESULTS 

The material combinations, the external loads, and the 
geometrical conf i gurati ons used in the numerical examples are 
summarized in Table 1 (see Figure 1 for the notation). 


Table 1. The cases considered as numerical examples 



The materials roughly correspond to a metal-hard polymer combi- 
nation (say, aluminum-epoxy) (cases A, B, and D) and an elastic 
half plane (cases C, E, and F). The crack is along the line 
( 0 = 9 Q , a<r<b) and the distances c, d, and the crack length 2a Q 
are defined by (Figure 1) 

2a Q = b-a , c = (a+b)/2 , d = c cos0 Q . (34) 

As indicated by Table 1, the results are obtained for uniform 
tractions 

T 2ee (r > e o> = Pi (r) “ ' p o • T 2re (r>e o ) = 0 ’ or 
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P-|(r) = 0 , P 2 ( r ) = ~P 0 » ( a< r < b ) (35. a, b) 

applied to the crack surface. The results are given in Table. 2. 
The quantities k.! , ( i =1 ,2) ' shown in the table are the stress 
intensity factor ratios defined by 

k’-(Cj) = k. (Cj)/(p p ^ap * (i ,j = 1 ,2; c^a, c 2 =b) (36) 

where p Q /a^ is the stress intensity factor in a homogeneous 
infinite plane with a crack of length 2a Q . The cleavage angle 
P c given in the tabl e' was obtai ned from (32). All results are 

'fc 

obtained for the plane strain case ( i . e . , k • = 3-4v . , j = l ,2) . 

J vJ 

In cases A, B, and C the distance d from the interface is 
fixed as d = 2a and the crack angle 0 Q is varied (see Figure 1). 
The limiting cases 6 Q = 0 and 0 Q = 90° of this problem were given 
in [7] and [4] and were reproduced with the present computer 
program for verification. Note that if the crack is in the less 
stiff material .( i . e. , >E 2 , case A), generally there is a 

reduction, and if E -j < E 2 (case B and C) there is an increase in 
the stress intensity factors compared to the values for the 
homogeneous medium. 

In cases D and E the radial distance c of the crack center 

is fixed as c = 2a and again is varied. In case D where 

o ,o 

E 1 / E g = 22.22 for 0 q = tt/ 2 crack becomes an interface crack for 
which the closed form solution is given by (e.g., [2, 3]: 


See [7] for the comparison of the plane strain and the plane 
stress results. 
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f 2 ( t ) + if, (t) ° °° b (,t°) (t - ZiB)R(t) , { 1 1 1 < 1 ) , 


iT 2r6> 


e 0 =n/2 


= (o, 


lT 0 )[(t - 21 B)R( t) - 1] , 


R(t) - (fry) 


- 1/2 


b r = 


n, + Kl v 2 


■» ■■ w 


b 2 ~ y 2 + K ? y 1 


t = (2r - b - a)/ (b-a) , 

k, - ik« = lim (too q ~ ( r _ a ) 

1 1 r->b+Q > oy ^ r0 0 =tt/ 2 r a 

o 

= (,% - i't 0 )(.l - 2i 3) /(b-a)/2 , 
~( a o + iT 0 } = + 1 p 2 ^ r ^ 5 


( 1 1 i >1), 


y = b -| / b 2 , 


U-j + K 1 y 2 


(b-a )/2 


(37) 


v/ here it should be noted that for the uniform tractions considered 
in (37) the stress intensity factors at r=a and r=b are the same. 
It is seen that in this case the stress singularity is oscillating 
in character and the definition of the stress intensity factors 
kj and k 2 is slightly different. Hence, at 8 = ir/'2 one would not 
expect the stress intensity factors and the cleavage angles to be 
continuous functions of -0 . The quantity which is expected to be 
continuous in 8 Q is the strain energy release rate given by (33) 
for a homogeneous medium. For the interface crack this quantity 
is given by [14] 




v (y-| + K-,y 2 ) (y 2 + K 2 y 1 ) 

2 y-]y 2 [(l+K-j )y 2 + 0 +k 2 'J 


i ; 2 2 

7TTI7I (k l + k 2> 


_ 1 + Ko 2 2 

5- — - — — ( kf + k 9 ) 

2 y 2 a 21 1 2 


(38) 
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us, for the crack imbedded in medium 2 if we evaluate the 
rain energy release rate ratio as 


»(H.c 2 )p 0 ‘a 0 2 


( 9U — V = k ' 2 + k 1 2 
l 3aJo 1 2 » 


( 39 ) 


jlrom (38) and 

iS / 3 U > _ / 3 U ^ 

Ji lim (tt - ) (57-) 

*■ 6 Q + tt/ 2 3a o 2 3 0 

shoul d have 


o 12 


(40) 


G 0 +tt/2 


( k-j 2 + k 2 2 ) = 


T~ W., 
a 21 21 


k 2 + k 2 
2 K 1 k 2 

a 21 P 0 2a o 


(1 + 3) 


(41) 


case D, a^-j = 3 . 93086 , 6 = 0.13420 and the limit becomes 

*, 0.51 796. which is given in Table 2. If one plots W 2 vs. 9 ( 
p may be S6en that t,iere 1 s in fact a smooth transition from 
laHe imbedded crack to the interface crack. 


|M In case E as 6 Q tt/2 the crack approaches the traction-free 
liu r face ; hence, as shown in the Table, the stress intensity fac- 
ors tend to infinity. It should be pointed out that the analysis 
iven in this paper is for a crack imbedded into a homogeneous 
edium. in the half plane problem for 0 Q = tt/2 the crack disap- 
ears, there is a discontinuity in the solution, and hence the 
f a)ues given in the table for 0 Q = tt/2 simply indicate the trend. 

The results given for case F show the effect of the radial 
sstance c (Figure 1) on the stress intensity factors k^ and the 
Jleavage angle (f» c for a fixed value of the crack angle, 0 o = 4O°. 
r - seen that as c increases the results approach the values 
.or ; ospondi ng to the infinite plane. 
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Table 2. The stress intensity factors and the probable 
crack propagation angles for the crack orientations, the 
material combinations, and the loads shown in Table 1. 
(k- = ki/Po/a^, i = l,2, a 0 =(b-a)/2). 


e o 

0 

20° 

o 

O 

60° 

CO 

o 

o 

90° 

k -[ ( a ) 

0.9349 

0.9307 

0.9216 

0.9144 

0.9143 

0.9160 

k ■] ( b ) 

0.9617 

0.9572 

0.9457 

0.9318 

0.9206 

0.9160 

k^U) 

0 

0.0025 

0.00001 

-0.0071 

-0.0153 

-0.0188 

k‘(b) 

0 

0.0125 

0.0209 

0.0237 

0.0215 

0.0188 

+c (a) 

0 

-0.306 

-0.001 

0.889 

1.911 

2.342 

^ c (b) 

0 

-1.492 

-2.531 

-2.915 

-2.666 

-2.342 

k ’ (a) 

0 

0.0223 

0.0338 

0.0333 

0.0251 

0.0191 

k ; C b ) 

0 

0.0030 

0.0020 

-0.0042 

-0.0142 

-0.0191 

k£(a) 

0.9349 

0.9395 

0.9492 

0.9580 

0.9644 

0.9663 

k^b) 

0.9617 

0.9629 

0.9655 

0.9671 

0.9672 

0.9663 

•» c (a) 

-70.53 

-70.08 

-69.85 

-69.87 

-70.03 

-70.15 

0 c ( b ) 

-70.53 

-70.47 

-70.49 

-70.61 

-70.81 

-70.91 

kj (a) 

1.0780 

1 . 0929 

1.1165 

1 .1389 

1.1459 

1 .1420 

kj(b) 

1 .0464 

1.0571 

1.0796 

1 .1091 

■ 

1 .1344 

1.1420 

k 2 ( a ) 

0 

-0.0098 

-0.0077 

-0.0062 

0.0242 

0.0321 

_ 

k 2 ( ) 

0 

-0.0220 

-0.0377 

-0.0433 

-0.0382 

-0.0321 

< M a > 

0 

1.024 

0.793 

-0.623 

-2.417 

-3.210 

<M b) 

0 

2.384 

3.985 

4.462 

3.852 

3.210 
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Table 2 (cont.) 


Case 

e o 

0 

20° 

o 

O 

■=3- 

60° 

o 

o 

00 

90° 

_ . . ■_ _U 

. 

k^ (a) 

0 

-0.0386 

-0.0586 

-0.0566 

-0.0399 

-0.0291 

■ 

k’ (b) 

0 

-0.0082 

-0.0092 

-0.0002 

0.0184 

0.0291 


k 2 ( a ) 

1.0780 

1 .0749 

1.0590 

1.0468 

1 .0400 

1.0390 

8 2 

k 2 ( b ) 

1.0464 

1.0465 

1.0418 

1 .0385 

1.0376 

1 .0390 


V (a) 

-70.53 

-71.22 

-71.59 

-71.56 

-71 .26 

-71.06 


<f c (b) 

-70.53 

-70.68 

-70.70 

-70.53 

-70.19 

-70.00 

■ . 

■ 

■ 

k-j (a) 

1 .0913 

1.1049 

1.1349 

1.1604 

1.1677 

1.1621 

' 

k ' (b) 
1 

1.0539 

T . 0644 

1.0918 

1.1257 

1 .1543 

1.1621 

r 

kj(a) 

0 

-0.0114 

-0.0093 

0.0067 

0.0276 

0.0367 

c i 

■■■■ 

k 2 ( b ) 

0 

-0.0251 

-0.0431 

-0.0498 

-0.0439 

-0.0367- 

• 

♦c< a ) 

0 

1.184 

0.941 

-0.663 

• 

-2.703 

• 

-3.61 

. 


V b > 

0 

2.700 

4.512 

5.041 

4.339 

3.61 


k i (a) 

0 

' 

-0.0440 

-0.0669 

-0.0646 

-0.0453 

-0.0331 


k-| ( b ) 

0 

-0.0095 

-0.0109 

-0.0003 

0.0205 

0.0331 

r 

k^(a ) 

1 .0913 

1 .0842 

1.0686 

1 .0540 

1.0452 

1 .0440 

L 2 

k 2 ( b ) 

1 .0539 

1 .0523 

1.0484 

1 .0444 

1 .0426 

1 .0440 


^c< a > 

-70.53 

-71 .31 

-71.73 

-71.70 

-71.36 

-71.05 


4> c (b) 

-70.53 

-70.70 

-70.73 

-70.54 

-70.15 

-69.93 
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if Tabl e 2 (cont. ) 


00 

0 

ro 

o 

o 

O 

O 

60° 

0 

o 

00 

kj (a) 

0.9349 

0.9201 

0.8749 

0.8034 

0.7334 

kj (b) 

0.9617 

I . . 

0.9523 

0.9217 

0.8631 

0.7753 

k 2 ( a ) 

0 

0.0018 

-0.0093 

-0.0539 

-0.1297 

k g ( b ) 

0 

[ .. 

0.0139 

0.0319 

0.0586 

0.0960 

4> c (a) 

0 

-0.227 

1.221 

7.608 

18.96 

*c (b) 

0 

-1.670 

-3.956 

-7.701 

-13.72 

W(a) 

0.8740 

0.8465 

0.7655 

0.6484 

0.5546 

W(b) 

0.9248 

0.9070 

0.8505 

0.7483 

0.6104 

kj (a) 

0 

0.0259 

0.0608 

0.1124 

s 

0.1547 

k-j ( b ) 

0 

0.0029 

-0.0008 

1-0.0258 

i , 

-0.0919 

k 2 ( a ) 

0.9349 

0.9299 

0.9160 

0.8910 

0.8012 

k‘(b) 

0.9617 

0.9586 

0.9493 

0.9305 

0.8598 

9 c ( a > 

-70.53 

-70.00 

-69.27 

-68.14 

-66.89 

0 c ( b > 

-70.53 

-70.47 

-70.55 

-71.06 

-72.58 

W(a) 

0.8740 

0.8654 

0.8427 

0.8064 

0.6658 

W(b) 

0.9248 

0.9188 

0.901 2 

0.8664 

0.7476 

k-j (a ) 

1 .0913 

• 

1 .1222 

' 

1.2437 

1 .6218 

4.0850 

k-j ( b ) 

1 .0539 

1 .0725 

1.1468 

1.3871 

3.0851 

k 2 ( a ) 

0 

j 

-0.0121 

-0.0041 

0.1125 

1 .4986 

kj><b) 

: 

0 

1 

-0.0287 

-0.0762 

-0.2099 

-1.2715 

■? c (a) 

0 

1.235 

0.373 

-7.859 

-33.458 

M b > 

0 

3.063 

7.540 

16.494 

36.013 




















Tab! e 2 (cont. ) 


Case 0 Q 


k\ (a) 
k] (b) 
k‘(a) 
l<2 ( b) 


v< a > 

^(b) 


0 

20° 

40° 

60° 

80° 

90° 

0 

-0.0524 

-0.1357 

-0.3262 

-0.9288 

— CO 

0 

-0.0104 

-0.0151 

0.0056 

0.1336 

CO 

1.0913 

. 

1.0979 

1.1177 

1.1381 

1.2538 

->■ CO 

1 .0539 

I 

1.0589 

1.0767 

i 

1.1187 

1.2702 

*->■ oo 

-70.53 

-71.44 

-72.86 

-76.08 

-84.91 


-70.53 

-70.72 

-70.80 

-70.43 

-68.53 


1.1 

1,2 

1.4 

1.6 

2.0 

4.0 

2.5263 

1.9959 

1.5947 

i 

1.4150 

1.2437 

1 .0543 

1.4694 

1.3772 

i ■ ■ 

1 1.2764 

i . 

1.2168 

1 .1468 

, ■- 

1.0423 

0.2103 

0.1091 

i 

0.0400 

1 

0.0141 

1 - 

-0.0041 

-0.0073 

-0.2911 

-0.2308 

-0.1625 

-0.1221 

-0.0762 

-0.0171 

• 

-9.388 

-6.219 

-2.870 

-1 .144 

0.373 

0.796 

. 

20.92 

18.08 

14.08 

11.24 

7.540 

1.882 

-1.0299 

-0.6667 

-0.3819 

-0.2538 

-0.1357 

-0.0227 

-0.0716 

-0.0491 

-0.0300 

-0.0219 

- 

-0.0151 

-0.0066 

1 .6185 

1 .4100 

I 

1.2580 

1.1884 

1.1177 

1.0288 

1.2639 

^ 1.2054 

1.1461 

1.1134 

1.0767 

1 .0227 

-82.94 

-79.73 

• 

-76.41 

-74.65 

• 

-72.86 

-70.95 

-71.61 

-71.31 

l 

-71.03 

-70.91 

-70.80 

-70.65 


























The results given in Table 2 may be used to obtain the 
tress intensity factors in an arbitrarily loaded two-phase 
omposite medium with an arbitrarily oriented crack provided 
he medium is loaded sufficiently far away from the crack region 
o that in the perturbed problem the crack surface tractions 
p^r) and p 2 (r) can be approximated by uniform stresses. Thus 
if p-|(r) + ip 2 (r) = p-j + ip 2 = constant, we find 


k-j ( c m ) * - I k ij( c m )Pj > (i ,j ,m = 1 .2; c^a , c 2 =b) 

(42) 

where in k!. i=l and i=2 respectively correspond to the normal 

' J 

and the shear components of the stress intensity ratio given in 

I | } 

Table 2 as k^ and k^, and j=l and j=2 respectively refer to the 
external loads p-j^O, p 2 =0 and p-j=0, p^O. For example, if the 
medium is loaded parallel to the interface away from the crack 
region, the uniform stresses in the uncracked material will be"' 
related by 

■ r i r0 ( r >°) = T 2ro(r,°) = 0 > 

1 - v, 2 1 - v 2 

Tj T ie0^ r ’°) = e ~ 2 T 2 e o ^ r ’ 0 ^ •* 


" 1 0 0 ( r , 0 ) 


T 200 (r ’ O) = °2 ’ 


(43) 


n r e a ^ and a 2 are constant. In this case the stress intensity 
'actors may be expressed as 

k i ( c ) ? , 

TJF' = k l (c m> = k il (c m> cos 0 o + k 12 (c m )s1 '"V ose o ■ 

-2‘ a o 


r 
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k p ( c„ ) , i p , 

— — - k 2 (c ra> = k 21 !c m )cos 6 o + k 22 (c m )si " e o c ° s8 o 
°2 / *o 


(m=l ,2 ; c-| =a , c 2 = b) 


( 44 ) 


figures 2-9 show some of the results obtained from (44). Figures 


2 and 3 give the results for the material combination and the 


crack orientation corresponding to case A in Table 1. Similarly. 


(aside from the external loads which are given by (43) ), Figures 


4 and 5 correspond to the case B, Figures 6 and 7 correspond to 


the case C, and Figures 8 and 9 correspond to the case D. Note 


that the probable cleavage angles cf> c shown in Figures 3, 5, 7, 


md 9 are all negative and the direction of crack initiation is 


approximately perpendi cul ar to the direction of the external 


oad. The angle <J> = 70.53 shown in the figures corresponds to 

v» 


■the cleavage angle for an infinite plane containing a crack and 


subjected to a uniform shear at infinity parallel and perpendic- 


lar to the plane of the crack [13]. For a fixed external load. 


figures clearly show the effect of the crack orientation on 


the stress intensity factors, and hence, on the fracture resist- 


■• rr 'o. of the composite medium 
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re 3. Cleavage angles and the strain energy release rates 
bonded planes shown in Figure 2 ( E, /E ? = 22 . 22 , v-, = 0. 

5,d=2a=b-a= constant). 













Figure 3. Stress intensity factor ratios in' bonded half planes 
containing an arbitrarily oriented crack and loaded parallel to 
the interface (E q /E 0 = 22.22 , v, =0.3, v 9 = 0.35, c = 2a = b-a = 
constant) . 1 ' . 





